arXiv: 1501.03131 v2 [q-bio.PE] 19 Oct 2015 


Hazard function models to estimate mortality rates affecting fish 
populations with application to the sea mullet (Mugil cephalus) 
fishery on the Queensland coast (Australia) 


Marco Kienzle* 
October 20, 2015 


Summary 

Fisheries management agencies around the world collect age data for the purpose of assessing the status of 
natural resources in their jurisdiction. Estimates of mortality rates represent a key information to assess 
the sustainability of fish stocks exploitation. Contrary to medical research or manufacturing where survival 
analysis is routinely applied to estimate failure rates, survival analysis has seldom been applied in fisheries 
stock assessment despite similar purposes between these fields of applied statistics. In this paper, we developed 
hazard functions to model the dynamic of an exploited fish population. These functions were used to estimate 
all parameters necessary for stock assessment (including natural and fishing mortality rates as well as gear 
selectivity) by maximum likelihood using age data from a sample of catch. This novel application of survival 
analysis to fisheries stock assessment was tested by Monte Carlo simulations to assert that it provided un-biased 
estimations of relevant quantities. The method was applied to data from the Queensland (Australia) sea mullet 
(Mugil cephalus) commercial fishery collected between 2007 and 2014. It provided, for the first time, an estimate 
of natural mortality affecting this stock: 0.22 ± 0.08 year” 1 . 
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University of Queensland, School of Agriculture and Food Sciences, St. Lucia, QLD 4072, Australia 
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1 Introduction 


One purpose of stock assessment is to estimate mortality rates affecting fish stocks. This estimation problem 
is easier to solve for species that can be aged as opposed to those for which age can’t be determined, for ex¬ 
ample crustaceans. The reason is that mortality and longevity are inversely related, hence age is a measure of 
mortality. The central mortality model in fisherie s research relati ng catch to the number of fish belonging to 


a cohort through time was proposed by Baranov (|Quinn and Deriso 


199911 . Given recruitment and mortality 


rates, the proportions of i ndivi duals at age in the catch can be calculated and used in a multinomial likelihood 


(Fournier and Archibald 


1982T) . This meth od has become by far the most common likelihood to integrate age 


data into modern stock assessment models (jFrancis , 


2014; 


Maunder and Punt 


20131) . 


The deterministic exponential model in Baranov’s catch equation has a statistical counterpart in the form 
of the exponential probability distrib ution fun ction which first and second moments quantify the relationship 
between longevity and mortality rate (ICowanl 1199811 : the mean age of a cohort which abundance declines at a 
constant rate is the inverse of that rate. Adopting such a statistical view of the exponential decay of individuals 
belonging to a cohort allowed the development of a set of maximum likelihood functions to estimate parameters 
of importance when assessing stocks. The field of survival analysis in statistics has created both a conceptua l 


framework and refined methods to estimate mortality rates (IKlcinbaum and Klein 
which are widely applied in medical research and engineering. 


2005 


Cox and Oakes. 


198411 


Despite the common goal of estimatin g survival rates i n medica l a nd fis heries research, survival analysis has 
seldom been applied to stock assessment. IFerrandis and Hernandea (120071 1 proposed to use the Weibull distri¬ 
bution as the survivor function to model data from scientific fishing surveys. In this manuscript, we developed 
an alternative application of survival analysis to model data from samples of commercial catches using hazard 
functions derived from time series of fishing effort and a schedule for gear selectivity. Likelihood functions of age 
data were derived to estimate a constant natural mortality rate, catchability and age-specific gear selectivity. 
This manuscript starts with a simplistic example using constant natural and fishing mortality rates to introduce 
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fundamental concepts from survival analysis applied to fishery research, before moving to more sophisticated 
cases leading to an application to real data from the sea mullet fishery in Queensland (Australia). The proposed 
methods were tested by Monte Carlo, using simulated data sets to characterize some of their properties and 
assert their capacity to estimate population dynamic parameters of interest to stock assessment. Finally, an 
application to the mullet fishery case study provided specific estimates of natural mortality, catchability and 
selectivity. 


2 Materials and methods 

Each fish can be assigned an age by examining its otolith, which is found just below its brain. Fish otoliths 
deposit calcium carbonate through time, thus increasing in size each year of their life. Microscopic observation 
of otolith sections often reveal alternate opaque and translucent zones, which can be used to assign individual 
fish to a particular age group. 


Sampling programs in fisheries research centers around the world aim to collect a representative sample of 
fish each year to determine the distribution of age of any species of interest. In most cases, the data are binned 
into age-groups of width 1 year. For this reason, we split the lifespan of cohorts from their birth (t 6 [0;oo]) 
into n yearly intervals from ai = 0 to the maximum age of a n+ 1 years. While the theory presented here used 
that particular subdivision of time (t), unequal ones also applies. 


2.1 The likelihood for constant natural and fishing mortality rates 


The exponential decrease in abundance of individuals belonging to a single cohort d ue to constant natural (M) 
and f i shing (F) mortality ra tes was described from a survival analysis point of view (IFerrandis and Hernandez , 


2007 


Cox and Oakes. 


19841) using a constant hazard function of time (t) and parameters 6 


h(t; 9 ) =M + F 


(1) 
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The probability density function (pdf) describing survival from natural and fishing mortality is 

f(t ; 6) = (M + F ) e ~( M+F)t = M x e ~ (M+F) * + F x e ~ (M+F) * (2) 

=/i(t;0) =h(t-,0) 

Since age data belonging to individuals dying from natural causes (note that contrary to human, fish’s largest 
cause of natural mortality is to be eaten by another fish) are generally not available to fisheries scientists, we 
used only the component of the pdf that relates to fishing mortality (f2(t;0)). This component of f(t;0) 
integrates over the entire range of t to 


pt — OO pt — OO 

/ / 2 (t;0) dt= Fx e~ {M+F)t dt 

Jt =o Jt =0 

pt — OO pt= oo 

= / fit] 6) dt — / M x e ~ {M+F)t dt 

Jt—o Jt =o 


/t=0 
= 1 - 

= 1 - 


M x e ~ {M+F)t dt 


t=o 

M 


M + F 

Hence, the pdf of catch at age data was obtained by normalizing / 2 (f; 0) 


(3) 


s(*;®) = 


M 


M+F 

M + F 


/ 2 (*; 0 ) 


F 

= mo) 


F x e" (M+F)t 


(4) 


The likelihood (Edwards 


19921) of a sample of fish caught in the fishery (Si) was written as 




n rt— a i 


f(t\9) dt) 


Si 


( 5 ) 


U P 


Si 


where Pi is the probability of dying in the interval [aj;aj+i]. 
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The logarithm of the likelihood was 


t—a. 


iog(£) = s i lo g( / /(*; d ) dt ) 

i= 1 Jt = a i 

n rt=GL %-\.1 

= 53 Si log( / (Af + F) e - (M+ ' F)t di) 

»=i ^ t= “‘ 

n 

= J2 s i log(e" (M+F)xai _e-( M+F ) x “ i+1 ) 


( 6 ) 


The log-likelihood can accommodate a last age-g roup made of all observations above a certain age in the 


sample (referred to as a -|-group) as follow (IPawitan . 


201,11 


log(£) = 53 Si log(e _(M+ ' F)xai - e -( M + F ) Xfli + 1 ) + l 0 g(e _(M+F)xOn ) 


(7) 


i=l 


This development illustrated an application of survival analysis to estimate mortality rates aff ectin g a cohort 
of fish by maximum-likelihood using a sample of catch at age. This method was implemented in R (IR Core Team . 


201311 in the package Survival Analysis for Fisheries Research (SAFR). 


Natural and fishing mortality cannot be disentangled with catch data only but the next section will show 
that the provision of effort data allowed to estimate both catchability (q) and natural mortality. 


2.2 Estimating catchability and natural mortality 

In this section, we assumed that a time series of effort (Ei) associated with a sample of catch at age (Si) was 
available to the researcher. And the assumption that fishing mortality varied according to fishing effort through 
constant catchability (q) held: F(t) = q E(t). In this situation, the hazard function was written as 


h(t,0) = M + q E(t) 


( 8 ) 


And the pdf 
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(. M + qE(t )) e - Mt ~lfo E W dt 

M x e - M ‘-9/o ®(‘) dt + q £(^) x £ (‘) dt 

s --.-' '-:-.-' 

=/i(t:6») =/ 2 (t;0) 


As in the previous section, we had 


(9) 


pt=c 

Jt =0 


/ 2 (t;0) dt = 1 - 


pt=c 

Jt =0 


M x e - Mt - q ti B V> dt dt 


( 10 ) 


But we did not know an analytic solution to the integral since the function E(t) was not specified. Never¬ 
theless, given effort in every interval (J t * ~^' +1 E(t) dt = Ei = f*Zo' +1 E(t) dt — J*~ 0 “' E{t) dt,Mi £ [l;n]), we 
could calculate the value of f *~ (3 °° f 2 (t; G) dt. 


pt — OO n 

/ / 2 (t;0) = 1-X)[ 

J t—0 i 


M 


M + q E, 


e -Mt-qJ*E(t) dty=<H +1 

J t=ai 




M 


M + q Ei 


,-M di-qf^ E(t) dt -M a i+1 -qf“' +1 E{t) dt\ 


( 11 ) 


_ ^ ^ (r~ M a i~ ( lfo i E G) dt _ P -M o i+ i-ij/ 0 a,+1 E(t) di^ 


' M + q Ei 
i=1 y 


In practice, // =0 °° / 2 (t; 0) is bound between 0 and 1. It takes a specific value depending on the values of M, g 
and Ei. Naming this constant value K , we could write the pdf of catch at age given effort data are available as 


g(t-,6) = ±Mt,0) (12) 

And the log-likelihood: 

n pt—ai -|-i 

log(£) = S) log( / g(t\ G) dt) (13) 

i =i “ ,t=a * 

Accounting for age-specific gear selectivity (s(f)) effects on fishing mortality (F(t) = q s(t) E(t)) was 
included in a similar way into the likelihood using constant value for selectivity at age. In practice, it is difficult 
to estimate n additional selectivity parameters using only the data from a single cohort but processing several 
cohorts at the same time assuming separability of fishing mortality rendered estimation of catchability, natural 
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mortality and selectivity possible. 


2.3 Estimates from catch at age matrix using fishing mortality separability 

This section describes an application of survival analysis to matrices of catch at age, developed for the purpose 
of estimating catchability (q), selectivity at age (s(t)) and constant natural mortality (M). The matrix (Sij) 
containing a sample of fishes aged to belong to a particular age-group j in year i contains n+p—1 cohorts. These 
cohorts were indexed by convention using k (k S [1, n+p — 1]) and an increasing number r /- (1 < < min(ra,p)) 

identifying incrementally each age-group (see Appendix p. [19] for more information). Each matrix Sij has two 
cohorts with only 1 age-group representing them. 


The derivation for a single cohort were the same as those presented in the previous section, here reproduced 
with indexations relative to a single cohort and accounting for selectivity 


9k(t-,0) = 


q s(t ) E(t) x e ~ Mt ~ q Jo s (*) E W dt 


\^ r k q £fc ,l E k j ( -M a k i-q L k ’ 1 s(t ) E(t) dt _ P —M a k t -q L k ' l+1 s(t ) E(t) dt\ 
2^1—1 M+q s k i E k i V ) 


(14) 


The likelihood function of a catch at age matrix was built using each pdf specific to each cohort (gk(t: 6)): 


n+p -1 r k ft—a k 1+1 

c = n nc/ ’ 9k(t-,e)dty Skd 

k -1 1=1 9t— a k,i 


The expression above is equivalent to 


£ = u< 


(15) 


(16) 


+3 


where the P t ,j are the probabilities of observing a fish of a given age j in year i given by the hazard model. In 
this likelihood, the P t ,j s um to 1 along the cohor t inste ad of summing to 1 for each year as described for the 


multinomial likelihood in 


Fournier and Arch ibald! ( 19821) . 


2.4 Monte Carlo simulations 

The method presented in the previous section to estimate mortality and selectivity from a matrix containing 
a sample of number at age were tested with simulated data sets to characterize their performance. Variable 
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number of cohorts (n + p — 1 = 25, 35 or 45); maximum age (p = 8, 12 or 16 years) and sample size of 
age measurement in each year varying from 125 to 2000 increasing successively by a factor 2 were used. The 
simulated data sets were created by generating an age-structure population using random recruitment for each 
cohort, random constant natural mortality, random catchability and random fishing effort in each year (Tab.[T|). 
A catch at age matrix was calculated using a logistic gear selectivity with 2 parameters: 


1 



1 + exp(a 

- 0 x ai) 

Variable type 

Distribution 

Parameters 

recruitment 

uniform 

min=10 b , max=l(y 

natural mortality 

uniform 

min=0.1, max=0.8 

catchability 

uniform 

min=10~ 4 , max=10 -3 

fishing effort 

uniform 

min=10 3 , max=5 x 10 3 

gear selectivity a 

uniform 

min=8, max=12 

gear selectivity 0 

uniform 

min=l, max=3 


(17) 


Table 1: Distribution and range of value taken by different type of random variable in simulations. 

Several sampling strategies were implemented to assess how it affected mortality estimates. To test estimators 
derived from survival analysis, one would like to draw randomly from the probability distribution. This is 
obviously impossible in the real world because field biologists never have in front of them a entire cohort to 
chose from. Nevertheless, we implemented a sampling strategy (sampling strategy 1) that randomly selected 
from the entire simulated catch at age matrix (C % j) as a benchmark. In the real world, samples can be drawn 
by accessing only a single year-class of every cohort every year, so the second strategy implemented was to 
simulate a random selection of a fixed number of sample ( N ) each year (sampling strategy 2). Finally, the third 
strategy investigated (sampling strategy 2 with weighting) was to apply a weighting by the estimated total 
catch at age ( Cij ) to the sample of number at age in the sample {Sij): 


Ci,j = Pij ®Ci® v(j) 


(18) 


where Pij is the proportion at age (see Appendix p. USD, c i is a column vector containing the total number 
of fish caught in each year i and v(j) is a row vector of l’s. A weighted sample (S'* J ) was obtained using the 
fraction of total catch sampled 


Stj = Cij x 


Yhi,j &i,j 

Ei Q 


(19) 







Note that &i,j = Eq; Stj- 


hj 


Comparisons with the multinomial likelihood proposed by I Fournier and Archibald (1198211 were made using 
differences in negative log-likelihood between that method and the hazard function approach described in the 
present article. Simulated catch were used to calculate the proportion of individual at age, constraining them 
to sum to 1 in each year. This method to calculate proportions for the multinomial likelihood was regarded as 
the best case scenario because we expect any estimation algorithm based on the multinomial likelihood to, at 
best, match exactly the simulated catch at age. The logarithm of these proportions were then multiplied by 
th e sim ul ated age sam ple ( weig hted or not depending on the case) to calculate the log-likelihood as described 


Fournier and Archibald! (1198211 . This quantity was compared to that calculated using the survival analysis 


approach to determine which model best fitted the simulated data. This comparison ignored the number of 
parameters used in each model. The multinomial likelihood requires n + p— 1 more parameters to be estimated 
than the survival analysis because the former requires an estimate of recruitment for each cohort in order to 
calculate the proportion at age in the catch. 

2.5 A case study: Queensland’s sea mullet fishery 

The straddling sea mullet (Mugil cephalus) is caught along the east coast of Australia, with most landings 
occurring between 19°S (approx. Townsville) and 37°S (roughly the border between New South Wales and 
Victoria). The most noticeable feature of the biolo gy of this spec ies is a massive northward spawning migra¬ 
tion of the stocks along the coast during autumn (IKcstcvenl . 1195311 . Ta gging experimen ts revealed that 90% of 


tagged animals travelled less than 85 km during the migration season (jKesteven 


19531 1. Analyses of parasites 


concluded that the bulk of sea mulle t caug ht in Queensland fishery is based on local fish pop ulatio ns and not 


migrating from New South Wales (ILester et all 1200911 . Following recommendations from 


Bell et al 


( 200511 . an 


existing (1999-2004) scientific survey design was modified from 2007 onward to include both estuaries and ocean 
habitats in order to provide representative demographic statistics of the fish caught in Queensland. Samples 
were collected in both habitats (Tab. GJ. Age varied between 0 and 17 years. A 14+ age-group was created to 
combine the small number of observations in the older age-groups. These data were weighted by catch in each 
year and habitat to obtain a dataset representative of the entire catch in this fishery. 


Sea mullet are thought to spawn in oceanic waters adjacent to ocean beaches from May to August each 
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year. By convention, the birth date was assumed to be on July 1 st each year. Opaque zones are thought to 
be deposited on the otolith margin during spring through early summer (September to December). Biologists 
have come to the conclusion that the first identifiable opaque zone is formed 14 to 18 months after birth, and 


all subsequent opaque zones are then formed at a yearly schedule (ISmith and Deguara , 


2003 ). Each fish in the 


sample was assigned an age-group based on opaque zone counts and the amount of translucent material at the 
margin of otolith. Age-group 0-1 comprised fish up to 18 months old (ai = 18 months) while all subsequent 
age-groups spanned 12 months (02 = 30 months, 03 = 42 months, etc ...). 


Three hazard function models were fitted to the data: a first model assumed a constant natural mortality 
across age-groups and throughout the period covered by the data, a common catchability and gear selectivity 
in estuaries and ocean (model 1, Tab. [3]); the second model assumed that catchability differed between estuaries 
and ocean; and the third model assumed that both catchability and gear selectivity differed between the two 
habitats. The models were compared usi ng Akaike_ Information Criteria (AIC) to determine which was most 


supported by the data (jBurnham and An ders on, 


20021 ) . 


3 Results 

3.1 Method tests using simulated data 

Weighting the numbers of sampled fish each year by total catch (sampling strategy 2 - weighted sample) per¬ 
formed as well as the benchmark sampling strategy 1 (Fig. Q] and Fig.[2]). By contrast, estimations using a fixed 
number of fish each year were biased suggesting that weighting by catch is necessary in practical applications 
of the survival analysis approach. 

Weighting of age-data samples considerably reduced the variability of natural mortality estimates (Fig |T]). 
Increasing the number of samples reduced uncertainty associated with natural mortality estimates too. 

Estimates of catchability were much more consistent across the range of values tested (3-10 x 10~ 4 ) for all 
methods (Fig. 0. The bias of the unweighted approach (strategy 2) was often similar to that of the weighted 
one (strategy 2 - weighted sample). But the uncertainty associated with the former approach was much larger 
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sampling.type 

Ej Strategy 1 
E-3 Strategy 2 

E-3 Strategy 2 - weighted sample 


Figure 1: Comparison between simulated natural mortality (x-axis) and estimated (y-axis) using a random 
sample of the matrix of catch (strategy 1); a random sample from each year separately (strategy 2) and the 
same sample weighted by yearly total catch (strategy 2 - weighted samples). Each panel correspond to an 
increasing number of samples per year varying from 125 to 2000. 
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Figure 2: Comparison between simulated catchability (x-axis) and estimated (y-axis) using a random sample of 
the matrix of catch (strategy 1); a random sample from each year separately (strategy 2) and the same sample 
weighted by yearly total catch (strategy 2 - weighted samples). Each panel correspond to an increasing number 
of samples per year varying from 125 to 2000. 


12 














than the latter. For both strategy 1 and strategy 2 with weighting, the benefit of increasing sampling size were 
very noticeable up to a 1000 fish aged but less so beyond that. 

The comparison between the likelihood function from survival analysis and the multinomial likelihood (Fig. [3]) 
showed that, apart sampling strategy 2 which provided biased estimates, the approach using survival analysis 
provided in the majority of cases smaller negative log-likelihood values than the multinomial likelihood. The 
substantial advantage given the multinomial likelihood in this comparison played an important role at low 
sampling intensity where the assumption that proportion at age was known perfectly artificially improved its 
performance in most difficult situations. This artificial advantage faded away as the simulated sample sizes were 
increased resulting in the survival analysis approach outperforming the multinomial likelihood. 


3.2 Mortality estimates for sea mullet 

Sea mullet data showed larger catch per unit of effort in the ocean than in estuaries (Tab. [2]). Of all three models 
compared with AIC, the model that assumed catchability varied between habitats and selectivity was the same 
in both habitats (model 2) was best supported by the data (Tab. El . This model estimated catchability in the 
ocean to be 16 times larger than in estuaries (Tab. QJ. Natural mortality for sea mullet was estimated to be 
equal to 0.219 ± 0.082 year -1 . Estimates of gear selectivity suggested it increased up to the fifth age-group, 
beyond which fishes were fully selected by the fishing gear. 


4 Discussion 

This application of survival analysis to fisheries research provided an effective approach to develop maximum 
likelihood estimators of natural and fishing mortality rates, and gear selectivity, from age data. Monte Carlo 
simulations showed that it provided unbiased estimates of natural mortality and catchability over a wide range 
of simulated values. 


The compari son between the negativ e 


mial likelihood (Fournier and Archibald 


log- likelihood from the survival analysis approach with the multino- 


1982T) suggested that the former offered a better model to represent 
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Figure 3: Difference between the negative log-likelihood (—log(£)) from survival analysis (SA) and multinomial 
(ML) as a function of the number of sample per year. Each panel represents a particular sampling strategy. 
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the data. This comparison was made using the best possible outcome for the multinomial likelihood because 
it used the simulated proportions of individuals at age in place of the probabilities to compute the likelihood. 
Arguably, a substantial advantage was given to the multinomial likelihood over the survival analysis in this com¬ 
parison because no one would reasonably expect any estimation method to systematically provide exactly the 
proportion at age in the catch using a sample of the data. Therefore, the present comparison really focused on 
which probabilities to use in the likelihood function, whether they should sum to 1 in each year along age-groups 
or along cohorts. Despite the strong advantage given to the multinomial likelihood, the results suggested that 
simulated data according to Baranov’s catch equation were fundamentally better fitted by a statistical method 
that modelled the exponential decay of individuals along cohorts rather than by one that assumed the data 
followed a multinomial probability distribution specific to each year. 

Weighting of the sample provided unbiased estimates of natural mortality and catchability. Mortality esti¬ 
mates, in particular fishing mortality, depended on the magnitude of catch. The unrealistic sampling strategy 
which assumed that all catch data would be in front of the experimenter at once for sampling, accounted 
automatically for variation of catch and effort in each year because the abundance of each age-group in the 
catch determined the probability to choose at random an individual belonging to any age-group. In practical 
application of survival analysis to fishery research, weighting is necessary because one cannot know a priori the 
magnitude of catch in coming years. 

The Monte Carlo simulations used a logistic gear-selectivity to generate and fit the data although we would 
have preferred to generate data from a wide range of possible gear-selectivity functions or even using non- 
parametric procedures. Simulations showed that gear selectivity were the most difficult parameters to estimate. 
The sea mullet case study was in fact not fitted with a logistic curve but selectivity were estimated through 
a tedious process to search each proportion retained at age that best fitted the data as measured by the like¬ 
lihood. This process could not be automatized into the simulation testing framework to provide automatic 
identification of gear-selectivity. This aspect of the analysis was left out of the present manuscript for future 
work. Criticisms that this somewhat simplified the problem would be justified. But the current article was 
designed as an introduction to the application of survival analysis to fisheries catch at age data, not one that 
solves all problems at once. As such, the likelihood approach presented in this manuscript provides a method 
to identify the gear selectivity that best fit the data, just not an automatic one. 
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The model best supported by the mullet data set estimated natural mortality equal to 0.22 ± 0.08. This is 
the first estimate of natural mortality for mullet in Australia. Previou s to th is estimation, it was customary to 
use the natural mortality estimated by linear regression from Hwang (119821) for the mullet fishery in Taiwan 
(M=0.33 year -1 ) which fall within 2 S.D. of the estimate for the Queensland fishery. The model that fitted 
best the mullet data estimated catchability in the ocean to be 16 times larger than in estuaries. This is consis¬ 
tent with fishermen reporting very large catches from their ocean beach operations (up to 40 tonnes per haul) 
compared to working in estuaries. 


This likelihood method may well find its place naturally into integrated stock assessment ((Maunde r and P unt 


20131) as it provided an efficient method to deal with samples of age data. Applications of survival analysis to 
fishery data could be expanded further. A particular area of interest would be to derive recruitment estimates 
using the probabilities estimated by survival analysis and total catch from the fishery. 
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Table 2: Distribution of yearly samples (in rows) of sea mullet into age-groups of width 1 year (in columns) in the estuaries 
and ocean habitats; catch in tonnes and effort in number of days. 
















































Model 

P 

-log(£) 

AIC 

2 

17 

18795.3 

37624.6 

3 

31 

18787.7 

37637.4 

1 

16 

18817.7 

37667.4 


Table 3: Comparison between the log-likelihood value obtained for several hazard function models of the mullet 
data using different numbers of parameters (p). The models were ordered by increasing values of Akaike 
Information Criteria (AIC =—21og(£) + 2 p) from top to bottom. 


Parameters 

Estimates 

Qestuaries 

(4.527 ± 1.026) xl0 b 

Qocean 

(7.243 ± 1.142) xlO” 4 

M 

0.219 ± 0.082 

Sl 

0.002 ± 2 x 10" 6 

S2 

0.08 ± 0.011 

S3 

0.374 ± 0.016 

S 4 

0.771 ± 0.045 

^5 

0.993 ± 0.07 

se 

1.000 ± 0.056 

sr 

1.000 ± 0.021 

S8 

1.000 ± 0.111 

S9 

1.000 ± 0.195 

Sio 

0.905 ± 0.326 

$11 

0.904 ± 0.503 

Sl2 

0.905 ± 0.389 

Sl3 

0.905 ± 0.528 

Sl4 

1.000 ± 0.661 


Table 4: Maximum likelihood parameters estimates from model 2. 
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5 Appendix 


5.1 Definitions of some mathematical symbols 

This appendix contains definitions of some of the mathematical symbols used in previous sections 


• Sij: a matrix of dimensions n x p (i £ [l,n] and j € [l,p]) containing a number of fishes that were aged 
and found to belong to specific age-groups j in a particular year i. This matrix contains data belonging 
to n + p — 1 cohorts, which by convention were labeled using k varying from 1 on the top-right corner of 
the matrix to n + p — 1 on the bottom-left (Tab. [5]). 


1 


1 


n 


n + p — 1 


j -P 

..321 
. 3 2 

. 4 3 

k . 


Table 5: Convention used to associate each element of the catch at age matrix (C,;j) with particular cohort 
referred to as with the number given in this table. 


The number of data in S, j belonging to each cohort (rfc) varies from 1 to min(n,p) and was determined 
as follow: 


Tk 


i — j + p if k < min(n,p) 

min (n,p) if min (n,p) < k < ma x(n,p) 

j — i + n if k > max(n,p) 


( 20 ) 


Each element of the Sij matrix is uniquely identified using indices i and j ( 1 < i < n and 1 < j < p) or 
indices k and l ( l<k<n+p — 1 and 1 < l < ), so for example 


= ( 21 ) 


• pij: a matrix of dimensions n x p (i G [1, n\ and j € [l,p]) containing the proportion at age in the sample 
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(Sij). Rows of this matrix sum to 1. 


Pi ,3 = 




Yhn &i,j 


( 22 ) 


• F t _j a matrix of fishing mortality with dimension n x p (i S [l,n] and j £ [l,p]). This matrix was 
constructed as the outer product of year specific fishing mortality rates (q Ei) and selectivity at age (sj): 


Fi t j =q Ei® Sj 


(23) 
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